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SCIENTIFIC  BACKGROUND  AND  RESEARCH  GOALS 


We  do  not  currently  understand  what  determines  the  fractional  cloudiness  in  partly  cloudy 
boundary  layers,  how  it  is  influenced  by  cloud-top  entrainment  instability  (Randall,  1980),  or 
what  controls  the  transition  from  stratocumulus  to  cumulus  conditions.  These  questions  have  a 
direct  bearing  on  the  radiation  budget  at  both  the  top  the  atmosphere  and  the  sea  surface,  and  also 
on  the  entrainment  rate  (Randall,  1987). 

The  importance  of  boundary-layer  clouds  stems,  to  a  large  extent,  from  their  powerful 
influence  on  the  net  radiation  at  the  sea  surface.  They  block  incident  solar  radiation  (with 
albedoes  up  to  about  50%  )  and  at  the  same  time  emit  downward  strongly  in  the  infrared.  Their 
net  effect  on  the  surface  energy  budget  is  believed  to  be  a  strong  reduction  in  the  energy  absorbed 
by  the  ocean  (e.g.  Esbensen  and  Kushnir,  1981).  Despite  the  obvious  importance  of  this  effect  for 
the  thermal  structure  of  the  upper  ocean,  very  little  is  known  about  the  coupling  between  the 
ocean  mixed  layer  with  boundary  layer  clouds.  The  coupling  can  work  in  both  directions;  the  sea 
surface  temperature  strongly  influences  the  type  and  amount  of  boundary-layer  cloud,  while  large 
cloud  amounts  are  favored  where  the  sea  surface  temperature  is  low  (e.g.,  Hanson,  1990). 

Modeling  studies  of  boundary-layer  clouds  date  back  to  the  work  of  Lilly  (1968),  who 
proposed  a  highly  innovative  yet  simple  model  of  the  cloud-capped  atmospheric  boundary  layer. 
Although  his  model  allowed  only  fully  overcast  or  completely  cloud-free  boundary  layers,  he  did 
point  out  the  possibility  that  the  cloudiness  could  be  reduced  by  various  processes,  including  what 
has  come  to  be  called  "cloud-top  entrainment  instability"  (CTEI).  Lilly  also  recognized  that  the 
cloud  amount  tends  to  be  larger  over  cold  water,  and  smaller  over  warm  water. 

Kraus  and  Leslie  (1982)  combined  an  atmospheric  boundary-layer  model  similar  to  Lilly's 
with  a  simple  upper  ocean  model.  The  joint  evolution  of  the  lower  atmosphere  and  upper  ocean 
was  predicted  by  the  model.  The  cloudiness  predicted  by  the  atmospheric  model  modulated  the 


insolation  of  the  upper  ocean,  thus  regulating  the  sea  surface  temperature.  Kraus  and  Leslie 
pointed  out  that  stratus  clouds  tend  to  keep  the  sea  surface  cool  by  blocking  solar  radiation,  while 
at  the  same  time  cool  sea  surface  temperatures  are  favorable  for  the  formation  and  persistence  of 
stratus  clouds.  They  recognized  the  importance  of  increasing  sea  surface  temperatures  for 
reducing  the  cloudiness,  although  their  model  was  not  able  to  predict  fractional  cloud  amounts. 

Randall  (1980,  1987,  1989)  has  recently  developed  a  simple  atmospheric  boundary  layer 
model  cq)able  of  predicting  the  boundary-layer  cloud  amount.  This  opened  up  the  possibility  of 
extending  the  study  of  Kraus  and  Leslie,  by  constructing  a  coupled  ocean-atmosphere  model  in 
which  the  cloud  amount  varies  dynamically  according  to  the  predicted  sea  surface  temperature 
(and  other  factors).  Such  a  model  can  be  used  to  study  the  transition  from  stratus  to  shallow 
cumulus  clouds,  and  the  role  of  air-sea  interactions  in  determining  the  time  and  location  of  such 
transitions. 


To  provide  the  data  needed  to  analyze  these  processes,  the  Office  of  Naval  Research  is 
sponsoring  ASTEX,  the  Atlantic  Stratocumulus  Transition  Experiment,  which  will  be  conducted 
near  the  Azores  in  the  summer  of  1992.  The  data  that  ASTEX  collects  will  be  particularly  useful 
if  they  can  be  applied  to  test  suitable  theoretical  and  numerical  models  of  the  marine  boundary 
layer.  A  primary  goal  of  this  project  has  been  to  develop  models  that  can  be  tested  in  ASTEX.  A 
further  goal  has  been  to  participate  in  the  ASTEX  planning  process.  Our  objectives  are  to  produce 
an  extensive  set  of  theoretical  and  numerical  tesults,leading  to  better  physical  understanding  of 
the  cloudy  marine  boundary  layer  and  providing  a  theoretical  basis  for  the  planning  and  execution 
of  ASTEX  (the  Atlantic  Stratocumulus  Transition  Experiment).  In  particular,  we  have  worked; 


•  To  gain  an  improved  theoretical  understanding  of  and  predictive  capability  for  partly 
cloudy  bound^  layers,  and  to  test  these  ideas  against  data  acquired  in  the  field. 


To  investigate  the  role  of  air-sea  interactions  in  regulating  cloud  amount  in  the  marine 


boundary  layer. 
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SCIENTIFIC  RESULTS: 

Marine  Boundary  Layer  Model  Development 

Randall  et  al.  (1992)  have  developed  a  new  type  of  boundary-layer  niodel  that  combines 
second-order  closure  with  a  bulk  representation  of  the  vertical  structure.  The  boundary-layer 
depth  and  turbulence  kinetic  energy  (TKE)  are  prognosdcally  determined.  The  large  turbulent 
eddies  that  are  primarily  responsible  for  the  fluxes  are  modeled  as  convective  circulations,  with 
ascending  and  descending  branches.  The  interior  of  the  boundary  layer  is  bounded  above  by  a 
thin  entrainment  layer  and  below  by  a  thin  ventil^on  layer.  Conservative  variables  such  as  the 
equivalent  potential  temperature  have  quadratic  profiles  in  the  interior.  Convective  circulations 
occur,  with  rising  branches  occupying  fractional  area  s,  which  is  predicted  by  the  model.  The 
upper  ocean  is  represented  by  a  mixed  layer  whose  depth  can  be  either  fixed  or  variable, 
depending  on  the  objectives  of  the  numerical  experiment  being  conducted. 

The  new  theoretical  model  has  been  developed  and  used  to  derive  the  scientific  results 
described  below.  The  same  model  is  being  tested  as  a  numerical  prediction  tool.  We  have 
completed  extensive  tests  of  the  model  against  two  data  sources: 

•  observations  collected  by  tethered  balloon  on  San  Nicolas  Island  during  FIRE  1987; 
and 

•  large-eddy  simulations  performed  by  Chin-Hoh  Moeng  of  NCAR. 

We  have  also  used  both  data  sources  to  evaluate  certain  parameters  that  appear  in  the  model.  This 
work  has  been  accepted  for  publication  in  the  Journal  of  the  Atmospheric  Sciences.  In  addition, 
we  have  coupled  the  boundary-layer  model  with  simple  upper  ocean  and  sea  ice  models.  We  are 
currently  investigating  the  possiblity  of  multiple  equilibria  of  the  ocean  -  ice  -  atmosphere  system, 
made  possible  by  the  boundary-layer  clouds  which,  when  present,  strongly  affect  the  energy 
budget  of  the  sea  ice. 
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Fig.  1  shows  an  example  of  results  obtained  with  this  model,  and  a  comparison  with  large- 
eddy  simulations. 


Change  of  q  with  height 


qu  and  qd  (g/kg) 


Figure  1:  Moisture  profiles  in  updrafts  (qu)  and  downdrafts  (qd)  as  simulated  by  a  laige-eddy 
model  (dashed  lines)  and  as  pr^icted,  based  on  fluxe  profiles,  by  the  bulk  boundary 
layer  model  developed  here.  The  larger  values  are  for  the  moist  updrafts,  and  the 
smaller  values  are  for  the  dry  downdra^. 

The  model  has  successfully  predicted  the  vertical  profiles  of  the  updraft  and  downdraft  properties 
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(in  this,  case,  mixing  ratio)  except  that  the  downdrafts  very  close  to  the  surface  are  predicted  to  be 
drier  than  the  large-eddy  simulation  predicted. 

The  new  boundary  layer  model  makes  a  number  of  new  predictions  concerning  the 
convective  turbulence  of  cloudy  layers.  For  the  special  case  of  a  well-mixed  layer,  it  predicts  that 
the  fractional  area  covered  by  rising  motion  is  near  1/2,  and  that  dissipation  in  the  interior  of  the 
layer  is  weak.  When  the  dissipation  is  weak  and  the  fractional  area  covered  by  rising  motion  is 
small,  the  nwdel  gives  the  "compensating  subsidence  -  detrainment"  relationship  that  has  become 
familiar  in  cumulus  parameterization  theories.  When  the  dissipation  is  strong  and  the  fractional 
area  covered  by  rising  motion  is  near  1/2,  the  model  gives  downgradient  diffusion.  For  the 
shallow  cumulus  regime,  the  model  predicts  that  the  fractional  area  covered  by  rising  motion  is 
smaller  for  the  case  of  large-scale  rising  motion  than  for  large-scale  sinking  motion.  These  model 
predictions  are  consistent  with  a  variety  of  observed  balances  in  convective  layers.  The  model  is 
based  on  various  assumptions  which  have  been  tested  against  data.  A  particularly  interesting  new 
result  is  the  prediction  that  potential  temperature  should  increase  upward  in  a  convective  layer, 
while  mixing  ratio  should  decrease  upward.  Both  predictions  are  consistent  with  observations. 

These  results  provide  a  very  useful  dynamical  interpolation  between  the  "compensating 
subsidence"  and  "mixing  length"  regimes.  They  tie  together  ideas  that  were  previously  thought  to 
be  umelated.  This  unification  opens  the  door  to  new  predictive  tools  that  are  both  simpler  and 
more  powerful  that  those  in  current  use.  These  ideas  have  already  been  useful  in  the  scientific 
planning  of  ASTEX.  In  addition,  we  have  demonstrated  that  the  predictions  of  the  model  are 
consistent  with  both  observations  and  LES  results.  This  provides  a  basis  for  confidence  in  both 
the  simple  riKxlel  and  the  large-eddy  model. 

It  is  our  intention  to  run  the  model  in  the  field,  in  “real  time”  numerical  experiments 
during  ASTEX. 
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Investigation  of  physical  processes  in  the  marine  boundary  layer 

In  a  further  investigation  of  the  effects  of  various  physical  processes  on  the  turbulence  in 
the  cloud-topped  boundary  layer  (Moeng  et  al.,  1992),  we  analyzed  three  large-eddy  simulations 
of  idealized  stratus-topped  boundary  layers  to  study  four  physical  processes;  cloud-top  longwave 
radiative  cooling,  entrainment,  surface  heating,  and  latent  heating.  Within  highly  turbulent, 
convective  boundary  layers,  turbulent  circulations  can  be  well  characterized  by  the  updraft  and 
downdraft  branches,  and  turbulent  transports  of  heat  and  moisture  can  be  well  described  by  the 
differences  between  these  two  branches. 

We  conditionally  sampled  the  large-eddy  simulation  held,  and  investigated  the 
contribution  of  each  process  to  the  mean-held  differences  between  updrafts  and  downdrafts. 
Based  on  the  results  of  this  analysis,  we  process-partitioned  the  total  heat  hux  and  the  total 
moisture  hux,  which  are  linear  in  height  for  a  homogeneous,  quasi-steady  state.  The  heat  and 
moisture  huxes  obtained  by  combining  the  partitioning  huxes  agree  quite  well  with  those 
obtained  directly  from  the  large-eddy  simulations. 

Fig.2  shows  that  the  updraft-downdraft  model  succeeds  very  well  in  reproducing  the 
large-eddy  simulated  hux  prohles  for  a  two  different  cases,  including  the  clear  convective 
boundary  layer  and  a  nocturnal  cloud-topped  boundary  layer.  The  agreement  is  amazingly  good, 
and  is  equally  good  in  two  more  cases  that  are  not  shown  here. 
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Figure  2:  Flux  profiles  obtained  through  large-eddy  simulation  (solid  curves)  and  as  predictecT 
by  the  updraft-downdraft  model  (dashed  curves),  for  two  different  physical  situations: 
a)  a  noctural  marine  boundary  layer;  b)  a  clear  convective  boundary  layer. 

Coupled  ocean-atmosphere  model  development 

In  addition,  we  have  constructed  a  coupled  ocean-atmosphere  model  which  can  be  used  to 
study  the  interactions  between  boundary-layer  clouds  and  the  sea  surface  temperature.  The  model 
currently  is  one-dimensional,  and  has  a  simple  Richardson-number  dependent  mixing  in  both  the 
atmosphere  and  ocean.  Mean  vertical  motion  and  geostrophic  shear  can  be  prescribed  for  both  the 
atmosphere  and  ocean.  It  is  our  intention  to  add  sea  ice  as  an  option  in  the  model,  in  the  future. 

Our  first  tests  were  for  cloud-free  cases.  These  were  used  essentially  for  debugging  the 
mixing  and  vertical  advection  codes.  An  example  is  shown  in  Fig.  3.  Here  we  performed  two  one- 
day  simulations. 
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Figure  3:  Results  from  a  test  of  the  coupled  ocean-atmosphere  model.  In  the  left  panel,  solar 
radiation  is  assumed  to  be  absorbed  completely  in  the  top  layer  of  the  model.  In  the 
right  panel  it  is  allowed  to  penetrate. 

In  the  first,  solar  radiation  was  assumed  to  be  deposited  completely  in  the  top  layer  of  the  ocean 
model,  while  in  the  second  the  solar  radiation  was  allowed  to  penetrate  through  several  layers. 
The  results  show  greater  surface  warming  in  the  absence  of  solar  penetration,  and  greater 
warming  at  depth  when  penetration  occurs. 

We  are  currently  testing  the  model  with  a  “dry  cloud”  that  produces  radiative  cooling  but 
no  latent  heat  effects.  This  problem  was  first  studied  by  Lilly  (1968)  and  has  proven  itself  to  be  a 
valuable  testbed.  We  plan  to  continue  this  model  development  woric  by  incorporating  latent  heat 
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effects,  and  adding  a  sea  ice  parameterization. 
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